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without side-band addressing. We show through numerical calculations that only local modes can 
be excited to achieve entangling gates through shaping the laser beams, so that the complexity of 
the quantum gate does not increase with the size of the system. These results suggest a promising 
approach for realization of large scale fault-tolerant quantum computation in two dimensional traps 
architecture. 
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I. INTRODUCTION 

The trapped ion system is one of the most promising 
candidates for quantum computation In the system, 
each qubit is formed by the long lived internal atomic lev- 
els of each ion, and laser beams have been used to manip- 
ulate quantum states of ionsjl, Q- Different theoretical 
schemes have been proposed for quantum gate operations 
with the aid of the ion's collective motions 2, 3, 4, 5, G,'?! 
and many have been experimentally demonstrated with 
small numbers of qubits [1, [l^, [ll|, . 

To make a quantum computer have practical value, 
scalability of proposed architectures to many qubits is of 
central importance. However, scaling up to tens of qubits 
in ion-trapped system proves to be a great challenge and 
many efforts have been done in this area. The main dif- 
ficulty of the scalability is that the sideband addressing, 
which is believed to be necessary in the gate operations 
proposed in most schemes, becomes impossible as the ion 
number increase. There are two promising approaches in 
literature to overcome this scalar difficulty. One is to di- 
vide the trapped-ion quantum computer into operation 
region and memory region. Quantum gate operations are 
always implemented at the operation region where only a 
few ions form a small linear ion chains, while ions can be 
shuttled between the operation and memory regions [l^ . 
The other promising approach is to use appropriately 
shaped laser pulses to realize high-fidelity conditional 
quantum gates, where all phonon modes are considered 
in the gate operations even in a very large linear ion 
crystal 1J| . This approach may be significantly strength- 
ened in an improved scheme [15| where the tightly con- 
fined transverse modesQ are used to mediate quantum 
gate operations and an anharmonic trap potential is pro- 
posed to arrange the ions in a uniformly-spaced linear 
crystal. In this approach, the gate operation has a local 
character so that the complexity of the operations does 
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not depend on the chain size [6*1 . 

To date most studies are limited to system of the lin- 
ear trap with ions laser-cooled and confined in one di- 
mensional crystals, little attention has been payed to the 
two-dimensional crystals[ii,[i3,[i3. Today penning trap 
can confine a large number of ions (10'* — 10^) by a po- 
tential with approximate cylindrical symmetry [l9|, |20| . 
In penning trap, axial confinement is realized by a static 
electric field, while radial trap potential is induced by 
the rotations of the ions under an axial magnetic field. 
Under certain parameters the ions will be arrayed in a 
two-dimensional plane with separation of order of tens of 
microns, such that they are individually addressable by 
optical means. Although the rotation of the ion crystal 
may induce additional technical difficulty in the address- 
ing of the ions, quantum manipulation of such qubits 
is still possible. Particularly it has been proposed that 
the axial phonon modes, which are independent on the 
phonon modes in the planar direction since the anhar- 
monic effects of phonon modes are negligible, can be 
used in the conditional gate operations 16]. However, 
the phonon modes in a two-dimensional crystal are much 
more complicated than that of a linear ion crystal, and 
such complex mode spectrum makes it more difficult to 
resolve individual modes for quantum gate operations. 
Clearly the two-dimensional penning trap may provide 
an almost ideal system for scalable quantum computa- 
tion and quantum simulation if the aforementioned diffi- 
culty of sideband addressing can be overcome. Therefore, 
a practical scheme to achieve the conditional gate oper- 
ations without sideband addressing, which has not been 
touched before in such large scale two-dimensional traps, 
is highly needed. 

In this paper, we propose a feasible scheme to imple- 
ment high-fidelity conditional phase gates without side- 
band addressing in the scalable two-dimensional traps. 
We investigate in detail the configurations of the ions as 
well as the normal modes of the two-dimensional crystals. 
This scheme uses the ion's motion along the axial direc- 
tion. Although the effects of phonon modes in the planar 
direction can be neglected, the axial phonon mode spec- 
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trum is still too complicated to be resolved in the gate 
operations. We show that the local modes suggested in 
Ref. [1, 0, [3] can also be applied to mediate the interac- 
tion between ion and laser beams, so that the complexity 
of the quantum gate does not increase with the size of the 
system. Compared to the scheme proposed in Ref(l6l|. 
we can implement conditional gate faster by optimally 
controlling the pulse shape. Since the nearest-neighbor 
entangling quantum gates are sufficient for fault-tolerant 
quantum computing in two dimensional trap, in sharp 
contrast to one dimensional chain where the nearest- 
neighbor entangling quantum gates are insufficient [23| . 
the generalization to two dimensional trap proposed here 
will be very useful. 

The paper is organized as follows: In Sec. II we study 
the equilibrium positions of the ions trapped in the two 
dimensional penning trap. In Sec. Ill we present the nor- 
mal modes of the trapped ions and especially investigate 
the properties of the axial modes in detail. In Sec. IV 
we show how to implement the conditional phase gates 
with high fidelity using the axial modes. A summary of 
our results and a conclusion are presented in Sec. V. 



II. EQUILIBRIUM POSITIONS OF IONS IN A 
TWO-DIMENSIONAL HARMONIC TRAP 

The s yste m we have in mind is a two-dimensional trap 
0, [13, l2ll | with N ions confined in a harmonic exter- 
nal potential. Including their mutual Coulomb repulsion 
interaction the potential energy of ions is given by 
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FIG. 1: (Color online). The Configuration of ion crystals for 
the number of ions N — 127. The circles represent the equi- 
librium positions of ions obtained by numerical calculation. 
For comparison, the solid squares denote the configuration of 
an ideal triangular lattice. Conditional phase gates acting on 
10 pairs of ions connected by arrows are studied. 



equations given as 
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where r„ = {xn,yn,Zn) denotes the position of the 
nth ion, M is the mass of each ion, e is the electron 
c harge, ep is the permittivity of free space and r„i — 
a/ {xn - XiY + {y„ - yi)'^ + {z„ - ziY is the distance be- 
tween the n-th and l-ih. ions. Here w,, (?] — x, y, z) is the 
trap frequency which characterizes the strength of the 
trapping potential in the rj direction. In the paper we as- 
sume that and [3 = Wx/w,.. The parameter 
/3 is a measure of the trap anisotropy and it will uniquely 
determine the ionic crystal structure for given N [22] . 

Similarly to what they do in Reffl^l for the ID linear 
case, to determine the equilibrium positions of ions, we 
may expand the position of the nth ion around its equi- 
librium position Yn^ as r„(t) w Yn^ -f q„(t), where qn(t) 
is a small displacement. The equilibrium positions are 
determined by the conditions 
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and its dimensionless form u™ = Ym /i (m = 1, 2, . . . , N) 
with i = ^e2/47reoMw2 

can be described by a set of 3A^ 
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where 

denotes the dimensionless distance between ions n and 
m. 

Determined by the number of ions and the trapped fre- 
quencies Ct'a:.^ z, the ions may form a string, a 2D or 3D 
ionic crystals [22|. To individually address qubits by laser 
beams, we usually assume ions are confined in a ID or 2D 
crystals. The ID case has been extensively studied, but 
less work focussing on the 2D traps[l^[ll|. In this paper 
we pay attention to the latter. Assume that lo^ is strong 
enough such that the configuration of ions in a 2D plane is 
always stable (the detailed requirement will be addressed 
below), then the equilibrium positions are confined in 
X — y plane and — for an arbitrary ion. Under this 
condition, the equilibrium positions 'mx—y plane may be 
determined by solving Eq.(l2]) and Eq.(l3|). The solutions 
are not unique, for instance, ions for = 4 may locate at 
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FIG. 2: (a) Dimensionless minimum separation as a function 
of total number of ions A'^ in a 2D harmonic trap, (b) The 
trap frequency tJr/27r (MHz) required for the minimum of 
distance dmin = 5/im, 20^m as a function of A^. Only the 
data for N = 7, 19, 37, 61, 91, 127 are given. 



the vertexes of a triangle and its center point or locate at 
the vertexes of a square. However, both experiments and 
theories show that the stable configuration of large num- 
ber of ions forms a distortional hexagonal latticejiol. |23 | , 
while the inner region arranges in equilateral triangles. 
In Ref.^3 laser-cooled ^Be+ ions in a two-dimensionally 
extended lattice planes were directly observed, especially 
a single plane hexagonal structure was observed, which 
is the configuration we want to utilize. For concreteness, 
in this paper we study the case where finite ions form 
the symmetric triangular lattice, such as the ion num- 
bers are = 7, 19, 37, 61, 91, 127, • • • . The configuration 
for N = 127 is shown in Fig. 1. In an ideal triangular 
lattice, the equilibrium position r^"^ may be expressed as 
r^°^ = [{j + l/2)aa, \/3ao/2] with denoting the lattice 
constant, j and I being integers. As expected the equilib- 
rium positions of ions obtained by numerical calculation 
exhibit little deflection from the ideal configuration. As 
seen from Fig. 1, the lattice spacings are not constant. 
As the ions get farther from the center, the lattice spac- 
ings become larger; however, it is almost constant near 
the crystal center and we get a nice triangular lattice at 
the center that can be suitable for the implementation of 
quantum computation. 

To manipulate quantum states of ions by laser beams, 
addressing individually the ions by optical means are 
needed, thus the minimum interior spacing must be larger 
than the size of the focal spot of the laser beam which is 
about 4/iTO in the experiments. The distances between 
the ions depend on the strength of the confining poten- 
tials as well as the number of ions. Actually the mini- 
mum value drain of the dimensionless distance between 
two adjacent ions descends as the total number of the 
ions is increased. By numerical fit from iV = 7 to 217, 
we find it obeys the relation Umin ~ 1.995/A^°-^''^, which 
is a slower rate than 2.018/jV"°-^^^[2l| in the linear chain 
case. The typical scaled distances M„im in 2D trap with 
small ion number N are shown in Fig. 2a. The relation 



between inter particle spacing dmin = Umin^ and the fre- 
quency ujr is tOr = \J 4^gp^j f — i SO ujr has an upper limit. 
The examples with small N are shown in Fig. 2b. 



III. NORMAL MODES 

We now turn to investigate the normal modes that 
are useful in the construction of the quantum opera- 
tions in ion-trapped quantum computation. The nor- 
mal modes in a linear trap has been studied in detail 
in Ref.[24], here we generalized the approach in one- 
dimensional crystal to two-dimensional penning trap. In 
a realistic ion trap the ions will have some nonzero tem- 
perature and will move around their equilibrium posi- 
tions. However, if the ions are sufficiently cold, we may 
expand the potential V under the harmonic approxima- 
tion (i.e., the terms O(q^) can be neglected), and obtain 

Trnlo- Here the subscript 
denotes that the double partial derivative is evaluated 
at q„ = q,„ = 0. In this case, the Lagrangian for the 
small oscillations is given by 
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The anharmonic terms in V will induce small errors and 
can be neglected, which has been proved in Ref . . The 
double partial derivatives in Eq.([5|) may be calculated 
explicitly, then the Lagrangian for the resulting small 
oscillations is found as 
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FIG. 3: (a)The minimum /3 required to confine all ions in a 
2D plane as a function of total number of ions in the Penning 
trap. Only the data for N = 7,19,37,61,91,127 are given. 
The line represents the best-fit behavior of the trend, (b) The 
energy gap of the axial modes between the center of mass 
mode and its nearest mode as a function of (5. 
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In term of the normal phonon modes, the motional 
Hamiltonian iJp for 3A^ motional modes can be written 
as the standard form 
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FIG. 4: (a)Quantum gate in a penning trap: electromagnetic 
fields induce a state-dependent dipole force on two ions in 
a triangular lattice, (b) A (T^-dependent force is driven by 
electromagnetic fields with two frequencies separated by /i, 
These fields couple the qubit states to the excited states with 
different coupling strengths, producing a differential AC Stark 
shift that oscillates at ^. The two fields must have a non-zero 
wavevector difference along the axial direction [26i]. 



that the 3D ionic structures become oblate to the x — y 
plane, forming a distorted 2D hexagonal lattice. In the 
following we assume that the parameter (3 is larger then 
the critical values. 

The spectrum of the axial modes exhibits some typical 
features: the center of mass mode at ujz is the highest fre- 
quency axial mode, while the frequency splitting between 
the center of mass mode and the second-to highest mode 
is the maximum, and the splitting will become smaller 
when increasing /3, as shown in Fig. 3b. We will impose 
the axial modes to implement quantum conditional phase 
gate in the next section. 



IV. CONDITIONAL PHASE FLIP GATES 



which includes three possible polarization: axial (z) 
modes and two in-plane (x, y) modes with the eigen- 
frequency . The axial and in-plane modes are indepen- 
dent, and the axial modes are determined by the matrix 
A^^ . The eigen- frequencies Wz^fc = ^JJT^ujr and eigenvec- 
tors b^'^^ of the axial modes are obtained from diagonal- 

ization of the matrix A^^ with Y.n.K^j^n'' = M^.fcb^'''- 
To make ions stable in the plane, all the angular fre- 
quency of the axial phonon modes i^z,/c, then the eigen- 
values iJLz,k of the matrix A^^{f]) must be positive. There- 
fore, the equation fiz,k = can be considered as the 
critical equation to determine the stability of the two- 
dimensional configuration. After solving this equation by 
numerical calculation, the critical values of /3 for different 
ionic number are shown in Fig. 3a. The trend can be well 
described by a simple pow law (3 = ^1.073(iV - 2)0-55, 
approximately coinciding with the results that others ob- 
tained when researching the structural transition of the 
crystalline confined system, such as the empirical scaling 
laws of Ref. [1^ , and the theoretic estimation of Ref. [l^l . 
So the critical value of /3 is the phase transition point 



We now start to address quantum gate operations 
in this two dimensional trap, mediated by many axial 
phonon modes. To perform quantum gates, we use the 
ac-Stark shift from two propagating laser beams with a 
relative angle and detuning to apply spin-dependent force 
on the ions aloiig the z direction(see FIG. 4) as they did 
in experiments [l^, [2^ . Under the Lamb-Dicke condi- 
tion the Hamiltonian of the system in the rotating-wave 
approximation becomes fsl. [l^ 

N 

H = -J2 W5^(4e--'=* + afce— (10) 

n,k—l 

where crfj = |0)(0| — |1)(1| is the Pauli operator, F„(t) 
is the spin-dependent force acting on the nth ion, and 
— \/h/2MuJz,ki>n^ is the coupling constant between 
the nth ion and the fcth phonon mode, b^''"' are the eigen- 
vectors of the matrix A^^ defined in the previous section. 
The spin-dependent force has the form of sin(/^t) deriv- 
ing from the ac Stark shift from the Raman laser beams, 
where fi is the detuning between the Raman laser beams 
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pair 
number 


U0r/2n = 0.2 


ujr/2n = 1.0 


F fi n^'-' 


F fi fi^" 


1 


0.999 10.033 0.177 


0.994 10.076 0.443 


2 


0.999 10.033 0.178 


0.989 10.076 0.442 


3 


0.999 10.033 0.173 


0.990 10.045 0.356 


4 


0.998 10.033 0.363 


0.989 10.045 0.443 


5 


0.998 10.033 0.169 


0.985 10.076 0.577 


6 


0.998 10.033 0.369 


0.973 10.025 0.534 


7 


0.997 10.033 0.566 


0.978 10.025 0.489 


8 


0.997 10.033 0.492 


0.979 10.025 0.505 


9 


0.997 10.087 2.768 


0.940 10.076 1.090 


10 


0.995 10.040 4.040 


0.690 10.076 1.315 



TABLE I: The gate fidelity F for 10 pairs of ions labeled in 
Fig.l with the same lJz(IOMHz) but two different ujr- The 
units for the detuning /i and the maximum Rabi frequency 
fl"' are l.OMHz. 

and 17 is the two photon Rabi frequency. The evolution 
operator determined by Hamiltonian pop can be expUc- 
itly found as [13 

U{t) e'*^" '^"''^^'^"+*^'<" (11) 

where (/)„(r) = Efc[an(T)al - a^i*(T)afc] with a'^ir) = 
fi lo ^n{t)g^e'^'^'-''*dt representing the residual entan- 
glement between ion n and phonon mode fc, and 

0/„(t) = ^ Jo dt2 EkFiit2)9f:9'nFn{ti) smu;,^kit2 - 
ti)dti characterizes the effective qubit-qubit interaction 
between ions I and n. Note that the phase 0/„(t) 
is actually a kind of unconventional geometric phase 
factor[27, 28], which has been experimentally demon- 
strated to have advantage of high-fidelity [!,[§]. 

Eq. pT|) is actually a conditional phase flip (CPF) gate 
acting on ions n and I if 0„ (r) = (pi (r) = and (pni (r) = 
7r/4. To implement a CPF gate on arbitrary ions n and I, 
we make the spin dependent force to be nonzero only for 
these two ions with F„(t) = Fi{t) — F{t), and set the gate 
time T so that (pniT) = (t>i{T) = and (j),ii{T) = 7r/4. In 
principle, it is always possible to satisfy these constraints 
by designing a sufficiently complicated pulse shape for 
the force. 

Taking the method in Ref. [l^, we use a simple se- 
quence of laser pulses which take minimum experimental 
control. A continuous-wave laser beam is chopped into 
m equal-time segments, with a constant but controllable 
Rabi frequency Qp for the pth {p = 1, 2, • • • , m) segment. 
The spin dependent force F{t) then takes the form F{t) = 
ripSin^fit) for the time interval (p — l)T/m <t< pT/m. 
With a small number of control parameters flp, a'^ir) 
(then 4'i{T)) in the evolution operation may not be exact 
zero, but as long as they are small enough, we still can get 
a high-fidelity gate. The task is to find an optimal control 
laser beam parameters sequence fip and /i to get a small 
enough infidelity Fin {= 1 — F with F being the gate 
fidelity). The gate fidelity is defined as F = 
where |^/) is the ideal output state after implement the 
CPF gate on the initial state j^I'o) which chosen typi- 
cally as l^-o) = (|0), + ® (|0)j- + \l)j)/2, and the 




FIG. 5: The gate fidelity F and the maximum Rabi frequency 
as a function of the laser detuning /i for ions pair 1 with 
LOr/2-K = 0.2MHz, uJ2tx = lO.OMHz and gate time 50^is. 

density operator pr = rr™(C/(T)(|«'o)(*o| ® p'^)U{t)'^) 
with tracing over the motional states of all the ions. Here 
we have assume that the phonon modes are initially in 
thermal state with an effective temperature T . 

Interestingly, we find that a very simple sequence of 
laser pulses ( such as m = 5) can ensure a good gate 
fidelity even for the two-dimensional trap. For concrete- 
ness, we calculate numerically the CPF gates between 
10 pairs of ions labeled in Fig. 1, with the same gate 
time r = 50/is and Wz/27r = lO.OMHz but two different 
LOr- The results are shown in Table I. In Table I, we list 
the optimal fidelity F, the corresponding detuning p and 
the maximum Rabi frequency f2*^ for each gate, where 
f2*^ — Max{f2p} can be used to characterize the maxi- 
mum laser power required in the gate operations. In the 
numerical calculation, we first calculate the fidelity and 
the required Rabi frequency for a CPF gate as a function 
of the detuning. As an example, the results for ion pair 1 
with ujr — 0.2MIIz is shown in Fig. 5, and for any other 
gate operations the relation between the fidelity and the 
detuning is similar, so fidelity has a local maximum near 
the detuning p — ujz- The data in Table I are derived 
from this approach. 

In Table I, the dimensionless distance between the 
pairs of ions is increasing in sequence. We observe sev- 
eral main features: (i) For both values of ujr^ it is clear 
that the gate fidelities decrease with the increasing of 
the distance between the pair of ions. Especially for 
bJr — l.OMHz, the fidelity decreases relatively quickly 
with the distance between ions, for example, the fidelity 
is about 0.994 for the 1st pair of ions, while it becomes 
0.690 for the 10th pair of ions. In addition, the laser pow- 
ers required for the CPF gates generally increase with 
the distance between the pair of ions, (ii) To understand 
the influence of minimum ionic separation to the gate 
fidelity, all gates are performed with two different trap 
frequency uJr- One frequency is set as ajr/'^n = 0.2MHz, 
and the other LOyjl-K = l.OMHz, then the corresponding 
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ion index 



FIG. 6: (Color online). The relative response of the ions for 
the gate operation, (a) The points (blank squares) in the x — y 
plane represent the ions in the equilibrium positions and the z 
values (filled circles) are the maximal deviation q^^J"' from the 
equilibrium position of the ions during the gate operation, (b) 
The response of the ions located along the three lines drawn in 
Fig. 1. The \q]^^\ for the target ions have been normalized to 
1. The fast decay of the response as one moves away from the 
target ions shows that the gate operation involves vibration 
of only local ions. 

minimum separation are 19.15/^to and 6.57/xm, respec- 
tively. From the table we can see that under condition 
of the same gate time and gate realized in the trap 
with larger minimum interior spacing will gain a higher 
fidelity. The physics behind it will be addressed late. 
On the other hand, it is worth to pointing out that, for 
pairs with lager distance, higher fidelity can be gained 
with lower Wj. or longer gate time, for example, taking 
gate time r = 200/is, the fidelity for the ion pair 10 with 
LOr/'i.TT = l.OMHz will bc above 0.995. 

The physics behind such high-fidelity is that the gate 
has a local character [l^, i.e., the contribution to the 
CPF gate comes primarily from the spin-dependent os- 
cillations of the ions that are close to the target ions. 
To demonstrate the local properties of the quantum 
gates, we plot in Fig. 6 the response of each ion dur- 
ing the gate operation between the 9th ion pair with 
w^/Stt = 0.2MHz. The figure explicitly exhibits the ultra 
fast decay of the ions response from the target ions. So 
the noises from the ionic vibration come just from the 
nearby ions. We have shown in Table I that the fidelity 
would be higher for a smaller since the ion spacings 
are less in the larger ujr- We have calculated the response 
of all ions for different w^, and find that of the nearby 
of the target ions for Ur — l.OMHz is larger of about 



10% than that of uJr = 0.2MHz. So one can under- 
stand that the fidelity would be lower by the increasing 
of the frequency tOr- Attributing to this local character 
of the gate operations, two notable advantages appear 
for the quantum computation: (i) the complexity of a 
gate operation in the two-dimensional traps does not de- 
pend on the system size, in sharp contrast to the schemes 
based on side-band addressing, (ii) The gates on the ions 
in different regions in the two-dimensional traps can be 
performed in parallel. 

The generalization of local gate operations from the 
linear trap0, fisj to two dimensional trap may have ad- 
vantages in the fault-tolerant quantum computation, ft 
is well known that the set of quantum gates with local 
entangling quantum gates and single-bit gates are uni- 
versal for quantum computation. Furthermore, it would 
be useful for large-scale fault-tolerant quantum comput- 
ing if the entangling quantum gates can be realized for 
any pair of qubits. However, although the quantum gates 
for any pair of ions in a large scale two dimensional trap 
are possible, they are slow and costly in the sense that 
the laser powers should be larger and the control param- 
eters should be increased if the distance between ions are 
large. In this case, an important question is whether it is 
still possible to perform large-scale fault-tolerant quan- 
tum computing with quantum gates that perform only 
on the ions with the distance much smaller than the size 
of the trap. The recent studiesjl^ have shown that in a 
two-dimensional lattice, even with only nearest-neighbor 
quantum gates, the error threshold can be still close to 
a percent level, which is basically as good as the case 
with arbitrary nonlocal gates, while the nearest-neighbor 
quantum gates are not sufficient for fault-torrent quan- 
tum computing in one dimensional chain. So the gen- 
eralization of trapped ion quantum computation to two 
dimensional trap is highly deserved. 



V. CONCLUSIONS 

In summary. We have explicitly exhibited the config- 
uration of the two dimensional trap and the equilibrium 
positions of the ions trapped in small ion number by nu- 
merical simulation, and then examined the condition on 
which all ions are confined in a plane. After investigating 
the properties of the axial modes of the trapped ions, we 
have shown through explicit examples and calculations 
that it is feasible to perform high fidelity CPF gate in the 
penning traps using the axial modes, and the gate has a 
local character as in the linear chain case. The results 
suggest a realistic prospect for realization of large scale 
quantum computation in the penning traps architecture. 
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